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We present a class of exact solutions of Einstein's gravitational field equations describing spher- 
ically symmetric and static anisotropic stellar type configurations. The solutions are obtained by 
assuming a particular form of the anisotropy factor. The energy density and both radial and tan- 
gential pressures are finite and positive inside the anisotropic star. Numerical results show that the 
basic physical parameters (mass and radius) of the model can describe realistic astrophysical objects 
like neutron stars. 
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I. INTRODUCTION 



Since the pioneering work of Bowers and Liang [1] there is an extensive hterature devoted to the study of anisotropic 
I spherically symmetric static general relativistic configurations. The study of static anisotropic fluid spheres is impor- 
tant for relativistic astrophysics. The theoretical investigations of Ruderman [2] about more realistic stellar models 
show that the nuclear matter may be anisotropic at least in certain very high density ranges (p > lO^^g/cm^), where 
^ ■ the nuclear interactions must be treated relativistically. According to these views in such massive stellar objects 
[ the radial pressure may not be equal to the tangential one. No celestial body is composed of purely perfect fluid. 
Anisotropy in fluid pressure could be introduced by the existence of a solid core or by the presence of type 3^1 su- 
' J ^ perfluid [3] , different kinds of phase transitions [4] , pion condensation [5] or by other physical phenomena. On the 

■ scale of galaxies, Binney and Tremaine [6] have considered anisotropies in spherical galaxies, from a purely Newtonian 
' point of view. Other source of anisotropy, due to the effects of the slow rotation in a star, has been proposed recently 
, by Herrera and Santos [7]. 

■ The mixture of two gases (e.g., monatomic hydrogen, or ionized hydrogen and electrons) can formally be also 

qh, described as an anisotropic fluid [8]. More generally, when the fluid is composed of two fluids the total energy- 

I ■ momentum tensor is 
5h 



r*^ = (Pi + pi) U'U^ - Pig''' + (P2 + P2) WW'' - P2<?'^ (1) 
where C/iC/* = 1 and WtW = 1. By means of the transformations 



> 
• 1— I 

>< 

' U*' ^Wcosa+J ^'^^^'^ W sin a, W*' = - J ^ W sina + W cos a, (2) 

]] Pi+Pi ]l Pi+Pi 

the energy momentum tensor (1) can always be cast into the standard form for anisotropic fluids, 

^ip + p) V'V" - pg'" + {a-p) (3) 

where V = U*'/./U*HJ*, x' = W*' / y/^W*HV* , p ^ T.kV'V'' , p = Pi + P2 and is a compHcated function of the 
densities and pressures of the two fluids [9] . 

The starting point in the study of fluid spheres is represented by the interior Schwarzschild solution from which 
all problems involving spherical symmetry can be modeled. Bowers and Liang [1] have investigated the possible 
importance of locally anisotropic equations of state for relativistic fluid spheres by generalizing the equations of 
hydrostatic equilibrium to include the effects of local anisotropy. Their study shows that anisotropy may have non- 
negligible effects on such parameters as maximum equilibrium mass and surface redshift. Heintzmann and Hillebrandt 
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[10] studied fully relativistic, anisotropic neutron star models at high densities by means of several simple assumptions 
and have shown that for arbitrary large anisotropy there is no limiting mass for neutron stars, but the maximum 
mass of a neutron star still lies beyond 3 — AMq. Hillebrandt and Steinmetz [11] considered the problem of stability 
of fully relativistic anisotropic neutron star models. They derived the differential equation for radial pulsations and 
showed that there exists a static stability criterion similar to the one obtained for isotropic models. Anisotropic 
fluid sphere configurations have been analyzed, using various Ansatze, in [9] and [12]- [19]. For static spheres in 
which the tangential pressure differs from the radial one, Bondi [20] has studied the link between the surface value of 
the potential and the highest occurring ratio of the pressure tensor to the local density. Chan, Herrera and Santos 
[21] studied in detail the role played by the local pressure anisotropy in the onset of instabilities and they showed 
that small anisotropies might in principle drastically change the stability of the system. Herrera and Santos [7] have 
extended the Jeans instability criterion in Newtonian gravity to systems with anisotropic pressures. Recent reviews on 
isotropic and anisotropic fluid spheres can be found in [30]- [31]. There are very few interior solutions (both isotropic 
and anisotropic) of the gravitational field equations satisfying the required general physical conditions inside the star. 
From 127 published solutions analyzed in [30] only 16 satisfy all the conditions. 

In the present paper we consider a class of exact solutions of the gravitational field equations for an anisotropic 
fluid sphere, corresponding to a specific choice of the anisotropy parameter. The metric functions can be represented 
in a closed form in terms of elementary functions. In the isotropic limit we recover the interior solutions previously 
found first by Buchdahl [24] and then by Durgapal and Bannerji [25]. Hence our solution can be considered the 
generalization to the anisotropic case of these solutions. All the physical parameters like the energy density, pressure 
and metric tensor components arc regular inside the anisotropic star, with the speed of sound less than the speed of 
light. Therefore this solution can give a satisfactory description of realistic astrophysical compact objects like neutron 
stars. Some explicit numerical models of relativistic anisotropic stars, with a possible astrophysical relevance, are also 
presented. 

This paper is organized as follows. In Section 2 we present an exact class of solutions for an anisotropic fiuid sphere. 
In Section 3 we present neutron star models with possible astrophysical relevance. The results are summarized and 
discussed in Section 4. 



II. NON-SINGULAR MODELS FOR ANISOTROPIC STARS 

In standard coordinates = {t, r, 6, (j)), the general line element for a static spherically symmetric space-time takes 
the form 

= A^{r)de - V-^{r)dr'^ - {d9^ + sin^ Bd^^) . (4) 
Einstein's gravitational field equations are (where natural units SwG = c = 1 have been used throughout): 

Ri - IrS^ = Tt- (5) 

For an anisotropic spherically symmetric matter distribution the components of the energy-momentum tensor are 
of the form 

Tt = {p + P±) Uivf' - p±6^ + (pr - p±) x^X^ (6) 

where is the four- velocity Au' = (5q, is the unit spacelike vector in the radial direction x' = VV5\, p is the 
energy density, Pr is the pressure in the direction of x* (normal pressure) and p± is the pressure orthogonal to Xi 
(transversal pressure). We assume pr ^ p±- The case Pr = p± corresponds to the isotropic fluid sphere. A = p± — pr 
is a measure of the anisotropy and is called the anisotropy factor [22] . 

A term appears in the conservation equations T^.^ = 0, (where a semicolon ; denotes the covariant derivative 

with respect to the metric), representing a force that is due to the anisotropic nature of the fluid. This force is directed 
outward when p± > Pr and inward when p± < p,.. The existence of a repulsive force (in the case p± > Pr) allows the 
construction of more compact objects when using anisotropic fluid than when using isotropic fluid [23] . 

For the metric (4) the gravitational field equations (5) become 



1-V V 2A'V V-1 

P= — 2 ,Pr = —r— + — —, (7) 

r Ar 



A ' r J ' \ A rA ) \ 
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where ' = 



It is convenient to introduce the foUowing substitutions [24]: 

rn(r) 



V=l- 2x7], x = r\v{r) = m(r) = i (0 d^- 

^ Jo 



(9) 



m(r) represents the total mass content of the distribution within the fluid sphere of radius r. Hence, we can express 
Eq. (8) in the form 



0. 



(10) 



For any physically acceptable stellar models, we require the condition that the energy density is positive and finite 
at all points inside the fluid spheres. In order to have a monotonic decreasing energy density P = {v'^) inside 
the star we chose the function r] in the form 



7] = 



ao 



(11) 



2(1 + V)' 

where tp = coa; and ao, co are non- negative constants. We introduce a new variable A by means of the transformation 



A = 



(ao - Co) (1 + V) 1 



ao 



(l-Ao)(l + V). 



We also chose the anisotropy parameter as 



3coAo-0 



(2 + Ao)(^ + l)' 



(12) 



(13) 



where Ao = ^ — 2. In the following we assume Ao > 0, with Aq = corresponding to the isotropic limit. Hence Ao 
can be considered, in the present model, as a measure of the anisotropy of the pressure distribution inside the fluid 

sphere. At the center of the fluid sphere the anisotropy vanishes, A(0) = 0. For small values of r, near the center, 
A(r) is an increasing function of r, but, after reaching a maximum, the anisotropy decreases becoming negligible small 
at the vacuum boundary of the star. 

Therefore with this choice of A, Eq. (10) becomes a hypergeometric equation, 



, (fA IdA 3 , „ 

^^-lhjT2 +o:7V-7^ = 0• 
dA2 2 dA 4 



(14) 



In the isotropic case Aq = wc obtain = §• Consequently, in this limit we recover the results obtained by 
Buchdahl [24] and Durgapal and Bannerji [25]. 

On integration we obtain the general solution of Eq. (14) and the general solution of the gravitational field equations 
for a static anisotropic fluid sphere, with anisotropy parameter given by Eq. (13), in the following form, expressed in 
elementary functions: 



^ = {(1 + i' f^ + /3i [5 - 2Ao + 2 (1 - Ao) V] ^2 - V' + Ao (1 + V)} 



(15) 



V = l 



3i/> 



(2 + Ao)(l + V)' 



(16) 



P = 



3co (3 + ^) 
(2 + Ao)(l + V)' 



-,Pr 



3co(p2 - Pi) 

P3 



,P± =Pr 



3AoCoV 



(2 + Ao)(l + V) 



2 ' 



(17) 



dpr 



3co 



P3 dtp \p3 



(18) 
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dp± dPr SAgCo (1 — V') 



# # (2 + Ao) (1 + ^pf ' 
where ai and /3i are constants of integration and we denoted 



(19) 



F (V-) = V2-^ + Ao(l + V), (20) 



pi = yrT? [3(^-1) -2Ao(l + V)] i=^(^^), (21) 



P2 = /3i 



3 (2^^ + 1) (V - 2) - 4Ag (1 + V)^ + 2Ag (TV'' + SV' - 2) - Aq (ISV-^ - TV' - 5) , (22) 



P3 = ( Ao + 2) (1 + V^) { (1 + V)'/' + /3i [5 + 2V^ - 2 Ao (1 + V)] (V') } (V') • (23) 

In order to be physically meaningful, the interior solution for static fluid spheres of Einstein's gravitational field 
equations must satisfy some general physical requirements. The following conditions have been generally recognized 
to be crucial for anisotropic fluid spheres [31]: 

a) the density p and pressure p^. should be positive inside the star; 

b) the gradients ^ and should be negative; 

c) inside the static configuration the speed of sound should be less than the speed of light, i.e. < ^ < 1 and 

o<^<i; 

d) a physically reasonable energy- momentum tensor has to obey the conditions p>Pr + '^P± and p + Pr + 2p± > 0; 

e) the interior metric should be joined continuously with the exterior Schwarzschild metric, that is A'^{a) = 1 — 2u, 
where u = M/a, M is the mass of the sphere as measured by its external gravitational field and a is the boundary of 
the sphere; 

f ) the radial pressure Pr must vanish but the tangential pressure p± may not vanish at the boundary r = a of the 
sphere. However, the radial pressure is equal to the tangential pressure at the center of the fluid sphere. 

By matching Eq. (16) on the boundary of the anisotropic sphere we obtain 

'^'■■' = ^- (2 + a!h1 + c„„^) =^-^"--- P") 

For the isotropic case, that is for Aq = 0, it is easy to show that 4(f^e°a^) ~ Therefore the mass-radius ratios 
for the anisotropic and isotropic spheres are related in the present model by 

2 

Uanis — 7^ ' 7 Uiso- (25) 
+ Ao 

Hence, the constants ai, /3i and cq appearing in the solution can be evaluated from the boundary conditions. Thus, 
by denoting X = coa' we obtain 

V 2 4«iso 2UaTOs(2 + Ao) 

X = coa = - — = - — ^ — —^-^ (26) 

ai = (1 - 2uanis) {(1 + Xf^ + A [5 - 2Ao + 2 (1 - Aq) X] ^2 - X + Aq (1 + X)}"' , (27) 



Q ^ V(l + X) [2-X + Ao(l + X)] [3(X-l)-2Ao(l + X)] 

^ 3(2X + l)(X-2)-4A3(l + X)V2A2(7X2 + 5X-2)-Ao(16X2-7X-5)' 
Note that equation (16) in [25] has been amended as the Eq. (28) presented here. 

In order to flnd a general constraint for the anisotropy parameter Ao, we shall consider that the conditions po = 
p (0) > 0, po = p{0) > and po > 3po hold at the center of the fluid sphere. Subsequently the parameters /3i and Aq 
should be restricted to obey the following conditions: 

< L iu..., Ao) = 3 + 2Ao-A(4Ag-4Ao + 3)V2TA^ ^ ^ 
- ^ ' ' l + /?l(5-2Ao)^/2TA^ " ^ ' 

The general behavior of the function L {uiso, Aq) is represented in Fig. 1. 
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FIG. 1. Variation of the function L (uiso, Ao) against the anisotropy parameter Ao and Uiso- 



Generally, the condition (29) is satisfied for < Aq < 1 and Uiso < 0.3. 



III. ASTROPHYSICAL APPLICATIONS 



When the thermonuclear sources of energy in its interior are exhausted, a spherical star begins to collapse under the 
influence of gravitational interaction of its matter content. The mass energy continues to increase and the star ends 
up as a compact relativistic cosmic object such as neutron star, strange star or black hole. Important observational 
quantities for such objects are the surface redshift, the central redshift and the mass and radius of the star. 

For a relativistic anisotropic star described by the solution presented in the previous Section the surface redshift Zs 
is given by 

/ 4 

The surface redshift is decreasing with increasing Aq. Hence, at least in principle, the study of redshift of light 
emitted at the surface of compact objects can lead to the possibility of observational detection of anisotropies in the 
internal pressure distribution of relativistic stars. Using the density variation parameter — p{R)/p{0), Patel and 
Mehra [29] discussed numerical estimates of various physical parameters in their model and concluded that the surface 
redshift in the isotropic case is greater than the surface redshift in the anisotropic case. Hence, our results are very 
similar to that of [29]. 

The central redshift Zc is of the form 

_ -1/2 ■ 



a 



1 + /?! (5 - 2Ao) Vs + A;^ -1- (31) 



The variation of the central redshift of the neutron star against the anisotropy parameter is represented in Fig. 2. 
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FIG. 2. Variation of the central redshift 2c as a function of the anisotropy parameter Ao for a static anisotropic fluid sphere 
with Uiso = 0.27. 



Clearly, in our model the anisotropy introduced in the pressure gives rise to a decrease in the central redshift. 

Hence, as functions of the anisotropy, the central and surface rcdshifts have the same behavior. 

The stellar model presented here can be used to describe the interior structure of the realistic neutron star. Taking 
the surface density of the star as ps = 2 x lO^^g/cm^ and with the use of Eqs. (17) we obtain 



or 



PsRn — 



Rn = 18.891w|/' (9 - Su.soY'" (2 + Ao) km, 



3X (3 + X) 
(2 + Ao)(l + X) 



2 ' 



1/2 



-1/2 



(32) 



(33) 



where i?jv is the radius of the neutron star corresponding to a specific surface density. For the mass and anisotropy 
parameter Ajv we find 



25.5684C (9 - Suiso^ (2 + Aq)""'" Mg 



1/2 



-3/2 



(34) 



Ajv = 7.19004 x lO^^AoWjso (9 - Suiso) ^ dyne cm~'^. (35) 

With Ao = we recover the results for {Mn,Rn) given by Durgapal and Bannerji [25]. In Eqs. (33)-(35), for 
the sake of simplicity, we have expressed all the quantities in international units, instead of natural units, by means 
of the transformations ^ 87r§f^, p p(? and A;v ^ ^Ajv, where G = 6.6732 x IQ-^dyne cm? g'"^ and 
c = 2.997925 x W^cm s"". 

The variation of the anisotropy parameter Ajv of the neutron star as a function of the radius Rn and mass Mjv is 
represented in Fig. 3. 
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FIG. 3. Variation of the anisotropy parameter Ajv (in units of lO'^^dynecm ^) of the neutron star as a function of the radius 
Rf4{km) and mass Mm (in solar mass units) for Ao G [0,0.9] and Uigo G [0.0001,0.28] . 

For a particular choice of the equation of state at the center of the star, Zpro ~ 3pxo = Po E^nd with a vanishing 
anisotropy parameter, Aq = 0, we obtain the result Ujso — 0.2908526 for a static isotropic fluid sphere that can be 
compared with the value of Uiso given in [25]. 

The quantities (yjj^^ , {j^^ i {^^^ (^^) ^"'^^ represented against the anisotropy parameter 

Ao in Fig. 4. 
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FIG. 4. Variations of the radial and tangential speeds of sound dp/dp at the vacuum boundary and at the center of the 
anisotropic fluid sphere as a function of the anisotropy parameter Ao for Uiso = 0.24: {dpr/dp)^^j^ (dotted curve); {dp±/dp)^^j^ 
(full curve); {dpr/dp)^^^ (dashed-dotted curve); {dp±/dp)^^Q (dashed curve). 

The plots indicate that the necessary and sufficient criterion for the adiabatic speed of sound to be less than the 
speed of light is satisfied by our solution. However, Caporaso and Brecher [26] claimed that dp/ dp does not represent 
the signal speed. If therefore this speed exceeds the speed of light, this does not necessary mean that the fluid is 
non-causal. But this argument is quite controversial and not all authors accept it [27]. 

IV. DISCUSSIONS AND FINAL REMARKS 

Curvature is described by the tensor field i? -^ j. . It is well known that if one uses singular behavior of the components 

of this tensor or its derivatives as a criterion for singularities, one gets into trouble since the singular behavior of 
components could be due to singular behavior of the coordinates or tetrad basis rather than that of the curvature 
itself. To avoid this problem, one should examine the linear and quadratic scalars formed out of curvature, such as 
ri =R,r2= R^JR'' and rg = R^juR''''' ■ 

With the use of the gravitational field equations (7)- (8) and of the static line element (4) we obtain the following 
expressions for the linear and quadratic scalars of the curvature tensor, given in terms of the radial pressure, energy 
density, mass and anisotropy parameter: 



n = 3pr - P + 2A, 



(36) 



r2= P 



2A (A + 2pr) , 



(37) 



r3 = 2A + p + ^ 



+ 2 



2m 



2m 



/2m 



(38) 



When A = p = = 0, we recover the invariant RijkiR^-'''^ for the Schwarzschild line element, that is ra = 48m^/r^. 
The variations of ri, r2 and rs are represented in Fig. 5. 
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FIG. 5. Variation of the curvature scalar ri (solid curve) and of the linear and quadratic scalars of the curvature ten- 
sor r2 (dotted curve) and rs (dashed curve) of the static anisotropic fluid sphere against the anisotropy parameter Ao for 
Uiso = 0.2908526. 

The scalars r2 and rs are finite at the center of the fluid sphere and monotonically decrease as the anisotropy 
parameter Aq increases. 

It is generally held that the trace T of the energy- momentum tensor must be non- negative. It is also the case that 
this trace condition is everywhere fulfilled if it is fulfilled at the center of the star [28] . 

The purpose of the present paper is to present some exact models of static anisotropic fluid stars and to investigate 
their possible astrophysical relevance. We have extended and generalized to the anisotropic case the method of 
obtaining exact solutions for relativistic spheres of Buchdahl [24] and Durgapal and Bannerji [25]. All the solutions 
we have obtained arc non-singular inside the anisotropic sphere, with finite values of the density and pressure at the 
center of the star. Variations of the physical parameters mass, radius, rcdshift and adiabatic speed of sound against 
the anisotropy parameter have been presented graphically. Our model can be used to study the interior structure of 
the anisotropic relativistic objects because it satisfies all the physical conditions and requirements (a)-(f). 
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